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Using the freedom in SO(10) GUT's one can generalize the existing models 
without changing the mass spectrum of fermions to obtain a significant suppression 
of proton decay resulting from the dimension 5 operators with AB ^ 0. In some 
limiting cases, these operators can be made negligible compared to dimension 6 
operators resulting from the heavy gauge bosons exchange. 
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1 Introduction. 



Proton decay is the most generic qualitative prediction of Grand Unified Theories (GUT's, 
see [[J , |L2|] , [|13|] , ]19| ) . Unfortunately, quantitative predictions are subject to several well- 
known uncertainties. Not only depend they on the choice of the GUT, but, for a given 
theory, serious uncertainties reside in various ambiguities in the values of the relevant 
parameters. In supersymmetric (SUSY) GUT's, the dominant contributions to proton 
decay amplitude usually result from the dimension 5, AB ^ operators ( || , [|iT| , |20|] , [|22| ) . 
describing effectively the exchange of a colour triplet with mass of the order of Mqvt'- 




The Wilson coefficients of the effective, dimension 5, AB ^ operators can be expressed 
in terms of the GUT parameters. The contribution of the diagram depicted above is: 



C 



abed 



1 



Mn 



■y a bycd 



Since these coefficients are suppressed with only one power of the large mass scaleQ, we 
can neglect the impact on proton decay of dimension 6 operators, induced by heavy 
gauge bosons exchange?]. The Wilson coefficients of the operators of dimension 5 have 
to be renormalized from the unification scale to the electroweak scale. At the scale of 
SUSY breaking, which lies rather closely to the electroweak scale Mz, the sparticles are 
integrated out and the operators of dimension 5 give rise to the effective four-fermion 
operators of dimension 6. This is depicted on the diagram below: 




1 In the minimal SU(5) Mt ~ Mqut but numerous SO(10) models admit Mr ~ Mpi anc k 

2 If the heavy gauge boson exchange is the only source of proton decay, the proton lifetime t ~ 10 35 y; 



see 16 
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When the gauginos are much lighter than the sfermions, the Wilson coefficients C of 
dimension 6 operators can be approximated by: 

16tt 2 m z - 

q,l 

Once the passage to the effective theory containing four-fermion operators is made, the 
renormalization procedure goes on to the scale of proton mass ~ 1 GeV. At this scale 
the link between three-quark operators and the interactions of hadrons is made using 
chiral Lagrangian technique, which exploits the approximate SU(3) y x SU(3) A flavour 
symmetry of the hadronic world (||, [/J). In this manner, proton lifetime can eventually 
be calculated. 

The first step (i.e. calculating the Wilson coefficients of the dimension 5 operators 
with AB ^ in a specific SUSY GUT) is relatively unambiguous in the simplest SU(5) 
model, which was widely investigated and served as a laboratory for developing the re- 
maining tools. However, the latest experimental lower bound on the proton lifetime in 
the dominant channel ( ||21|| ): 

r (p -> K + i>) > 6.7 • 10 32 y at 90% CL (1.3) 



indicates that the minimal SU(5) is ruled out (cf. |T7j for an updated discussion). There- 
fore, recently the attention has been focused on SO(10). Unlike the minimal SU(5), in 
which the Wilson coefficients of the d = 5, AB ^ operators are determined by the 
Yukawa couplings and Mt up to two phases^, SO(10) models are less constrained and 
thus less predicitive, but they give some hope to increase the prediction for proton life- 
time in a natural way. Recently, several results in SO(10) have been reported (H, |T%|). 



claiming maximal r v to be within one order of magnitude from the present experimental 
bounds. 

The purpose of this paper is to take a critical look at those predictions. We shall 
use some straightforward generalizations of the model presented in [|J to illustrate some 
features of SO(10) GUT's that allow to suppress proton decay induced by the operators 
of dimension 5. In Section 2 we discuss thoroughly the contributions to the Wilson 
coefficients of the operators of dimension 5 relevant for proton decay and compare them 
with the contributions to Yukawa couplings. In Section 3 we stress the importance of 
the proper renormalization procedure and point out some subtleties that led to errors in 
the published works. In Section 4 three strategies for increasing the prediction for the 
proton lifetime are proposed. In Section 5 we explain why in SO(10) GUT's one usually 
needs large mass scales of the order of the Planck mass, which does not happen in case 
of the minimal SU(5), and sketch another method for raising the proton lifetime. The 
conclusions are presented in Section 6. 

3 The impact of these phases on proton decay has been studied in 113]. 
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2 The link between fermion masses and proton decay. 



In this section we analyze a class of SO(10) GUT's which can provide a realistic mass 
spectrum for fermionsQ. In these theories the light Higgs doublets may reside in the vector 
representation 10h and in the spinor representation 16h, i.e.: 



H u — 2 10 



(u) 

H 



H d = cos 7 d2i a ; H + sin 7d 2^ H (2.1) 

We also assume that even at the GUT scale we are dealing with an effective theory 
resulting from Planck-scale physics which admits: 

• renormalizable interactions described by operators consisting of three superfields 
(Yukawa couplings) 

• non-renormalizable interactions described by operators consisting of four super- 
fieldsQ. They become effective Yukawa couplings, when one of the fields acquires a 
vacuum expectation value (vev) of the order of the unification scale. 

We also assume an additional global Z 2 discrete symmetry, under which 16 a are odd for 
a = 1, 2, 3 and all the remaining fields are even (which is called matter parity). 

We shall use the notation in which the symbol [...]#[...] signifies that the representa- 
tions in the left and the right square brackets are contracted to form a representation R 
(or R + R, if necessary), respectively. In this notation all the effective Yukawa couplings 
which are consistent with the assumed symmetries read: 



16 a 16 6 10 H 


(2.2 


[16 a 16 6 ] 10 [(45 H )10 H ] 


(2.3 


16 a 16 6 ] 120 [(45 H )10 H ] 


(2.4 


[16 (16 H )] 10 [16 6 16 H ] 


(2.5 


16 a (16 H }] 12 o[166l6 H ] 


(2.6 


[16 a (16 H )] 1 [16 6 16 H ] 


(2.7 


[16 a (16 H )] 45 [16 6 16 h ] 


(2.8 


16 a (16 H >] 21 o[16 6 f6 H ] 


(2.9 



4 See jj] for an example of a theory of fermion masses. 
Operators with more vevs have also been considered in the literature, but in a different context (Cf. 
p8[). We make a heuristic assumption that the effective couplings containing more vevs of the order of 
the unification scale are suppressed by higher powers of the quantity ^j GUT . 
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The symmetry group SO(10) can be broken to the Standard Model gauge group Gsm = 
511(3)^ x SU(2) L x U(l) y by the vevs of the adjoint and spinor Higgs fields: 

(45 H ) ~ B-L 

(16 H ) ~ v c (2.10) 

(the conjugate spinor representation 16h must also acquire a vev in SUSY theories). 
Then the coupling (|2.4j ) introduces mass differences between the down-type quarks and 
leptons, the coupling fl2.5|) and (|2.6|) provide a non-trivial Kobayashi-Maskawa matrix and 
the couplings ( p. 7| ) ( p79| ) generate Majorana masses for right-handed neutrinos when both 
the spinor representations 16h acquire vevs. 

Now we give the explicit formulae for the Yukawa couplings and the couplings to 
colour triplets resulting from the operators (|2.2|)-( [2T9"|) . We denote by T R (T R ) a triplet 
(antitriplet) coming from the representation R. After the GUT symmetry is broken, the 
part of the superpotential resulting from these operators reads^]: 

Wy = Y^ b U c aa U^H u + Y^ b & aa D^H d + Yf b E^ d + 

+y a l u° ^6^io e <*/3 7 + y^u^bifi^e^ + 

(u:E b - by b ) f 10a + {u:E b - by b ) f 16a + 

+M^A (2.11) 

The matrices Y u , Y D , Y E are the Yukawa couplings of the up- type quarks, down- type 
quarks and leptons, respectively. The Majorana mass matrix of the right-handed neutrinos 
is denoted by M R . The matrices and (AfW an d X®) are the Yukawa-like cou- 
plings of the left-handed (right-handed) matter fields and colour triplets and antitriplets, 
respectively. If we denote by the contributions resulting from eq. (|2.10| +n), all the 

6 We use Greek letters for colour indices, lower Latin letters i-r for weak isospin indices, lower Latin 
letters a-e and upper Latin letters for flavour indices. 
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abovementioned matrices can be expressed in the following way: 



Y ab 


(!) , 
= %b)+ a 


1 ab 


( m 

= -(<$> 


I ab 


- -(a {1) 




(!) i 
= %b) + a 




(7) , 

= a) (n + a 

(ab) 


n ] 


(1) 

~~ a (ab) a 




( 4 ) i 
= a (ab)+ a 


1,(1) 

ab 


(1) 

= %b)- a 


y(2) 
X ab 


( 4 ) i 
(ab) 


v(l) 
X ab 


(!) i 
= %b) + a 


ab 


(7) , 
= a (ab) + a 


Ml 


~ 1 a (ab) + 



(3) 
[ab] 



A3) 



.(4) 



,(5) 



(2) 
(ab) 

(8) 
(ab) 

(2) 
(ab) 

(5) 
(ab) 

(2) 
(ab) 

(5) 
(ab) 

(2) 
(ab) 

(8) 
(ab) 

5 

4 



- 2a 



(3) 
ab] 



+ 2a 



(3) 
ab] 



(7) 



(8) 



(ab) ^ "(ab) 



Mi 



(7 



(2.12) 
(2.13) 
(2.14) 
(2.15) 
(2.16) 
(2.17) 
(2.18) 
(2.19) 
(2.20) 
(2.21) 
(2.22) 

(2.23) 



where M v = |(16 H )|- 

Having integrated out heavy triplets, we obtain the following effective Lagrangian. 



f f i r^abcd/r\ . /~iabcd/r\ 

L — L.MSSM + t> Uabcd + t> L^abcd 



where: 



Cabcd — 
C^abcd = 



QaiQbjQlkLdl 

Ua a D c b ?U c n E< 



(2.24) 

(2.25) 
(2.26) 



We shall call (9 a bcd and O a bcd the LLLL and RRRR operators. 

We denote by M~ l the inverse of the mass matrix for the colour triplets. The elements 
of this matrix are labelled by the relevant representations of SO(10) to which the triplets 
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belong. Using the equations fl2.12| )-( |2.23p we obtain: 



?r abcd _ ( (1) , (2) \ / (1) _ (2) _ „ (3) \ M -l 

lLj - [ a (ab) + a (ab) ) [%d) a {cd) Za [cd]J iW 1010 + 

+ «i) + <$)) («$) - «S) - 2ag>) M^e + 
+ («8)+«8)) (^) + «!S))^o + 

+(«a-«S))( a S) +o S))^ i io+ 

+ (<$) + <$)) + a (2)) M w \ 6 (2.28) 
Factor 2 results from the following identity: 

QaiQbj £ ij £ <*i3-y = 2 Ufa D ^ (2.29) 

In the doublet-triplet splitting mechanism presented in [|J, Mj^ 16 = in order to keep 
one pair of the weak doublets massless at the unification scale. The remaining Higgs 
fields, including the colour triplets, have masses of the order of Mqut then. 

2.1 The freedom of a general model. 



The coupling ( p.3|) does not contribute to the Yukawa couplings but affects proton decay. 
Indeed, when the gauge symmetry SO(10) is broken to Gsm, the vector representation 
10h splits into two weak doublets which carry no B — L charge and two colour triplets with 
B — L = ±|. Because (45h) ~ B — L is evaluated directly on this vector representation, 
only the couplings to colour triplets can arise from 



The coupling (277) contributes only to the Majorana masses of the right-handed neu- 
trinos and has no impact on proton decay. This can be seen as follows. The vev (16h) is 
invariant under Gsm and it is contracted with some matter field to give a SO(10) singlet, 
which is, of course, a Gsm singlet. Thus, this matter field must be the only Gsm singlet in 
the matter multiplet: u c . Because both the right-handed neutrinos and the colour triplets 
decouple near the unification scale, the coupling fl2.7| ) cannot contribute to the dimension 



5 operators with AB ^ 0. Since the combinations of the couplings couplings (2.7)-( 279|) 
that enter the Majorana masses of the right-handed neutrinos and those enetring the 
amplitudes of proton decay are linearly independent, there is no generic link between the 
Majorana masses of the right-handed neutrinos and proton decay. 
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2.2 Babu-Pati-Wilczek model. 



In [|J, a quantitative GUT of the type described above is presented. A realistic spectrum 
of fermion masses can be obtained for = 1, = —0.1096, = 10~ 5 , — = 
c&zs = —0.0954, c&23 +«23' ) = —0.0413, +012'' = 4.14- 10~ 3 and the remaining entries in 
the matrices . . . ,a^ 5 ' equal 0. Some assumptions about the spectrum of the theory 



at the Planck scale (e.g. the couplings (|2.3| ) and ( |2.4j ) can only arise due to an exchange 
of 16 + 16 states with masses of the Planck scale, the matrices a^ 6 \ . . . , a*- 8 -* are chosen so 
that the contributions to the Yukawa couplings of neutrinos are equal to that to proton 
decay) remove all the freedom discussed above. This arbitrary choice makes the model 
maximally predictive, but the predictions for the proton lifetime are inconsistent with the 
experimental bound ( |1.3| ). Therefore, we shall treat this model a starting point for some 
generalizations which allow to increase the predictions for the proton lifetime. 



3 The importance of the proper renormalization. 

In the minimal SU(5) model, Wilson coefficients of the LLLL operators, which result from 



integrating out heavy colour triplets T, can be written as flT5f : 

1 



C 



((/)(■(/ 



ab -\-icd 



-y aa y 



where: 



y{M GUT ) 
y(M GUT ) 



Y U (M GUT ) 

^GUT) 



Y D (M C 



Y E {M GUT ) 



(3.1) 



(3.2) 
(3.3) 



Though the running of y and y differs from the running of the Yukawa matrices and 
the equalities ( |3.2| ) and ( |3.3| ) do not hold except at the unification scale, they can help 
to connect the observables from the Standard Model with the Wilson coefficients relevant 
for proton decay. The following two recipes can be applied. 

Recipe 1. Within the MSSM the Yukawa couplings y c (Mz), yd(Mz), etc. can be 
expressed by m c (lGeV), m rf (lGeV), etc using the renormalization group (RG) analysis. 
Next, Yukawa couplings are evolved upwards to M G ut, where y and y are calculated 
by using the equations ( |3.2j ) and ( |3.3|) . At the GUT scale we use ( |3.1| ) to calculate 
Wilson coefficients, which are then renormalized downwards to Mz- If a,b,c,d = 1,2, 
the renormalization group equations (RGE's) for C abcd and C abcd can easily be solved 
by using the 1-loop solutions for the running gauge couplings. We obtain the following 
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enhancement factors (J 



.4 



abed 



A 'abed 



C abcd {M z ) 
C abcd (M GUT 

gj(Mz) 

g!(M GUT ) 

C abcd (M z ) 
C abcd (M GUT 

gj(Mz) 

g£(M t 



gl{M z 

gl(M ( 



GUT j 



g\{M z 
gKM ( 



GUT j 



(3.4) 



GUT, 



g\{M z 
gI(M ( 



GUT, 



_2_ 

11 



(3.5) 



The RGE's for the Wilson coefficients of the operators involving the third generation 
of fermions must be solved numerically, since the running is affected by the large Yukawa 
couplings. At M z all SUSY particles are decoupled, and dimension 5 operators give rise 
to dimension 6 operators with coefficients C' abcd , which are renormalized downwards to 
1 GeV. This results in another enhancement factor: 



-4/ 



C 



abed i 



lGeV) 



C' abcd {M z ) 



1.35 



(3.6) 



Then the amplitudes of proton decay can eventually be calculated. 

Recipe 2. We can also use the equations (|3.2|) and (|3.3|) to define the quantity: 



Cf cd {t) 



Y a U b (t)Y c ^(t) 



Mn 



(3.7) 



which obeys: 



C: bcd (M GUT ) = C abcd (M GUT ) 



(3i 



This is not the Wilson coefficient relevant for proton decay, because it is renormalized in 
a different way. However, knowing the running of the Yukawa couplings and the Wilson 
coefficients, we can calculate the discrepancy between C® bcd (M z ) and C abcd (M z ). This is 
the suppression factor As introduced in ref. [jlll and calculated (for m t = 100 GeV) in 



ref. 0: 



.4c 



gj{M z 

gi(M ( 



GUTj 



g\{M z 
g!(M ( 



GUT) 



0.67 



(3.9) 



Then we can express Cf 3cd {t) in terms of the running masses. The coefficient of a prototype 
operator relevant for proton decay reads: 



C abcd {M z ) = A s Cf cd (M z ) = A 



abed i 



g\ 



ab-\rcd* 



M T M^ sin 2f3 H 



m ua (M z )m d c(M z )5 ab V c 



CKM 



(3.10) 
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since at the Mz scale all the Yukawa couplings can be expressed in terms of the running 
masses at that scale. After dressing with superparticles we can repeat this procedure 
to dimension 6 operators. Between Mz and 1 GeV the discrepancy between the exact 
solution of the RGE for true Wilson coefficient and the coefficients expressed by running 



masses results in another suppression factor introduced in refs. [11], |15|. Having solved 



the 1-loop RGE's for the strong coupling constant and the running masses, we obtain: 

0.43 (3.11) 



_ ( gj(m b ) yS (gl{m c )\-% ^ 3 2 (lGeV) \^ 
L \9i(Mz)J \gl{m b )) \ g!(m c ) J 



and the following relation: 



A~ L 6 = A L (3.12) 



It has been pointed out in ref. || (cf. a recent discussion in H) that this value differs from 
that given in the abovementioned papers. Solution of the 3-loop RGE's for the strong 
coupling constant and the running masses implies: 

A L = 0.28 (3.13) 

The relation ( |3.12j ) is also modified then. 

Note that the meaning of A s and Al is completely different from that of As and 
Al, respectively, since the former cover all the renormalization effects, not the discrepancy 
in running of two quantities. The Recipe 2 cannot be conveniently used in the SO(10) 
models, because, as we have demonstrated in Section 2, the Wilson coefficients relevant 
for proton decay cannot be parametrized with the Yukawa matrices in a simple manner. 
However, the coefficients A a s bcd and A s lbcd , which do not depend on the choice of GUT, 
can be calculated either in the minimal SU(5) or in SO(10). In the former case, they are 
functions of fermion masses at M^, tan/3 and the corresponding As factor^]. 



In particular, it is not legitimate to use the values of As and Al given in ref. ||T5[ to 
account for all the RG effects in SO(10) models. However, we have found that in some of 
the previous works on SO(10) models these renormalizations were computed incorrectly. 
For example, the authors of ref. M use the wrong value As = 0.67 and obtain the upper 
theoretical bound on the proton lifetime at least twice too large. 



4 How to increase the prediction for the proton life- 
time? 

In order to make definite predictions for the proton lifetime we assume that Majorana 

masses of the right-handed neutrinos are generated solely by the coupling ( |2.7|) and = 

7 In fact, As depends on the generation indices of the relevant Wilson coefficient and on the mass of 
the top quark. This has been analyzed in [JL5|. The value quoted here was obtained for m t = 100 GeV. 



9 



a® = 0, except for Section O , in which the possible link between the Majorana masses 
of the right-handed neutrinos and proton decay is investigated. We consider the following 
two SUSY spectra: 

• all sfermions have masses 1000 GeV except the lighter stop, whose mass is taken to 
be 400 GeV, the masses of charginos are 100 GeV and 500 GeV, gluinos have masses 
350 GeV 



• inspired by fl?], we also analyzed the case in which the first two generations of all 
sfermions have masses 10 TeV, the third generation has masses 1 TeV except the 
lighter stop with mass 400 GeV; we shall call this case the decoupling. 

We neglect the contributions of the neutralinos to proton decay. We also use the latest re- 
sults of the JLQCD collaboration for the hadronic matrix element \a\ = \(3\ = 0.015 GeV 3 
(p[). The decay channels are enumerated in the following way: 



No. 


decay channel 


1 


v - 


-> K+P 


2 


v - 


-> 7Y + U 


3 


v - 


-> K°fJL+ 


4 


v - 


+ 7TV+ 



4.1 Neutrino masses and proton decay. 

In order to study he impact of and a® on proton decay we assume that + 
is the fraction of the Majorana mass matrix of the right-handed neutrinos normalized to 



Mr, 



Mr, 



(4.1) 



We also assume = a®. In the limit of £ — > 0, the model becomes what is called the 
Strategy I with £ = 1 in the following section. The prediction for the proton lifetime as a 
function of ( for the two SUSY spectra is depicted on Figure |I[ The maximum is shifted 



towards ( < due to a destructive interference of the contributions resulting from ( |2.8| 
and (|2.9|) with the remaining ones. 



4.2 Strategy I. 

In [Q, in the model dubbed Case II, the coupling of the form 

[16 2 10h] 16+1 - 6 [(45 h }16 3 ] 



(4.2) 



10 



tr- 




37 
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34 
33 
32 
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Figure 1: The proton lifetime (in years) in the dominant channels as a function of (. Only 
the AB ^ operators of dimension 5 have been taken into account. The lower dotted 
line, the lower dashed line, the upper dotted line and the upper dashed line are the partial 
proton lifetimes in the channels 1,2,3 and 4, respectively. The solid line represents the 
current experimental lower bound on the proton lifetime in the decay channel 1. See 
Section [O] for more explanation. 
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is used to provide flavour-antisymmetric contributions to the Yukawa couplings which 
allow the down-type quarks and leptons of a given generation have different masses. The 
operator (|4.2|) is, of course, a combination of the operators fl2.3|) and (|2T4]) . This can also 
be obtained with the use of the coupling: 

-[16 3 10h] 16+1 - 6 [(45 h )16 2 ] (4.3) 

Then the operators ( [4.2|) and (|4.3| ) give the same contributions to a^ 3 \ whereas their 
contributions to a*- 2 ** have opposite signs. In general, both of them can be present in the 
theory. Then and are no longer equal but proportional: 

= £ a ® (4.4) 

The dependence of the proton lifetime on the parameter £ is depicted on Figures |2|. 

There are three special cases £ = —1,0,1, in which the sum of the operators (|2.3|) 
and (|2.4j) can be represented by a single SO(10) invariant operator. For £ = 1 there is 
an accidental cancellation of the various contributions to the LLLL operators. This is 
due to the fact that, numerically, y = Y U ^Y U and it is diagonalized simultaneously with 
the Yukawa couplings. The maximum is shifted to £ > 1 because of the interference with 
the RRRR operators. In the case of the decoupling the proton decay results mainly from 
the stop exchange in the SUSY loop and the coupling to the third generation is quite 
insensitive to changing £. 

4.3 Strategy II. 

It follows from (gg7p that if 



a® = -a« (4.5) 

the Wilson coefficients C abcd are driven to zero and only the operators consisting of right- 
handed superfields are relevant for proton decay. They are, however, strongly suppressed 
by the smallness of the Yukawa couplings of the first two generations of fermions. 

One could however argue that the assumption ([4.5|) is an extreme fine-tuning. There- 
fore, we shall consider a slightly more general scenario, in which: 

a< 2 > = £a« (4.6) 

The dependence of the proton lifetime on the parameter £ is depicted on Figure |[ 

For £ = — 1 there is a maximum of the lifetime in all the channels, because the 
contribution of the LLLL operators vanish and the proton decay is induced solely by the 
RRRR operators. The maximum at £ = 1 for the channels p — > K°fi + and p — > 7r°/x + 
results from vanishing the coefficients of the RRRR operators. 
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Figure 2: The proton lifetime (in years) in the dominant channels as a function of £ 
for Strategies I and II. See the caption of Figure (P and Sections fO| and fO| for the 
explanation. 
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The agreement with the experimental data can only be obtained in a narrow band 
around £ = — 1. Though this solution is slightly more natural in the case of decoupling, 
it still requires some fine-tuning of the parameters. 



4.4 Strategy III. 

This idea utilizes an additional adjoint representation 45r, which acquires a vev propor- 
tional to the right isospin of the Pati-Salam model ||19|| . Then the coupling (|2.2 ) can be 
replaced withf]: 

[16 a 10 H ] 16+1 - 6 [(45 R )16 o ] (4.7) 

Because the right isospins of the left-handed particles equal zero, this coupling does not 
contribute to C abcd . Moreover, the contribution to C abcd is antisymmetric in the first two 
indices, since the right isospin of the up-type antiquarks is minus that of the down-type 
antiquarks. 
We assume 

a& = £ fl < 3 > (4.8) 

The dependence of the proton lifetime on the parameter £ is depicted on Figures |||. 

The maximum of the proton lifetime in all the dominant channels at £ = results 
from vanishing the coefficients of the LLLL operators there and the prediction for the 
proton lifetime is consistent with the expreimental data. The case of £ = is equivalent 
to the absence of the coupling (|2.3|), so it satisfies the requirement of naturalness. 



5 Why does one usually need larger Mt in SO(10) 
models than in the minimal SU(5)? Strategy IV. 

In the minimal SU(5) there are two Yukawa couplings: one for the up-type quarks and a 
common coupling for the down-type quarks and leptons. The model predicts that at the 
unification scale the Yukawa couplings are equal to the relevant couplings to the colour 
triplets y and y (see equations ( |3.2p and (|3.3|)) and when the Yukawa couplings are 
diagonalized, y and y are determined up to two phases JT^|. After the renormalization 
to the electroweak scale and dressing with a SUSY loop, we have: 



(f**(M z ) ~ As ^* 2 m AM z )mAMz)J ab VS^ M ({U) 



CKM 

M T m\ M&ain2p H 



8 This pushes the value of V c b slightly out of the experimentally allowed region. However, we regard 
this strategy only an illustration of a general possibility of model building. 
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Thus, if the LLLL operators are dominant^: 

r P ~ ~ ( " T2 (5-2) 

and the proton lifetime is maximized for the models with small tan/?//. 

In many SO(10) models, the Yukawa couplings of the third generation of fermions 
result from the same coupling I63I63IOH, which makes them all equal at the unification 
scale (see e.g. fll8|). Since the top-bottom-tau unification requires large tan/5//, the eq. 
fl5.2|) yields that the proton lifetime is additionally suppressed in these models. 

There is also a class of models, in which the light doublet Hd is a mixture of states 
coming from the representations 10h and 16h as in (|2.1|) . In these models (e.g. |J, [|J): 

Vt tan/?// 

— = — 5.3 

Vb cos 7 d 

Thus tan (3 H can be small, provided that cos 7^ < 1. Assuming that the Wilson coefficients 
exhibit a hierarchy similar to that of the Yukawa couplings, we expect that the equation 
( p.ip becomes: 



C {Mz) ~ As M^\ M^sin2/?//cos 7d (5 " 4) 



(where all the quark masses are the running masses at Mz) and the smallness of cos7d 
additionally suppresses the proton lifetime: 

COS 2 7d /, r\ 

r P ~ ^2 ( 5 - 5 ) 

tan Ph + 



tan Pi 



This suggests another way of raising the predictions for the proton lifetime. In a model 
with small tan/?//, if we add the operators which make the Yukawa coupling of the top 
quark much larger than the other Yukawa couplings and do not contribute to proton 



decay, we could avoid the suppression (|5.5|) . Such a model has been proposed in |10| , but 
it reproduces the bad mass relations of the minimal SU(5). 

Despite this disadvantage, we analyzed a combination of the two models presented in 
||10|| . We assumed that the top quark is a mixture of the states coming from the spinor 
and adjoint representations: 

t = £«16 3 + a/1 - £ 2 «45 (5.6) 



9 In case of large tan/3//, the RRRR operators may become dominant. This does not spoil our ar- 
gument, since in that case r p ~ (tan/3//) -4 ( fl4|| ), so the proton lifetime is even shorter than in our 
estimation. 
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Then the large mass of the top quark is generated by the operator 45 2 together with the 
operators ( |2.5| )-(|2^9|)[ t0 |. Moreover, only ^y t instead of y t enters the relevant contributions 
to the proton decay amplitude. We also assumed that the large masses of the colour 
triplets are generated by the following part of the superpotential: 

W T = 6 1 16 H (45 R )r6 , H + &2l6 H (45 R }16i I + 6316^(1)16^ (5.7) 

where 16^ + 16' H are fields which neither couple to matter nor acquire vevs. Then the 
effective mass of the colour triplets: 

Mr = b -^f (5.8) 

can naturally be made larger than the unification scale Mqut- Therefore, we set Mt = 
2 ■ 10 18 GeV as in the other SO(10) models. The predictions for the proton lifetime as a 
function of £ are presented on Figure |3|. 

In all the cases, the predictions for the proton lifetime are consistent with the experi- 
mental data. 

6 Conclusions. 

We have analyzed a class of SUSY SO(10) GUT's which can provide fermions with realistic 
masses. We have shown there is a lot of freedom in the parameters of such theories which 
does not affect the predictions for the low energy phenomena except for proton decay. 
We have proposed four strategies to increase the predicted proton lifetime. Though our 
solutions are clearly lacking naturalness, they show that no definite quantitative results 
can be obtained for the proton decay induced by the d — 5, AB ^ operators. Therefore, 
only the predictions for proton decay induced by the heavy gauge bosons exchange, which 
are model independent, can be an unambiguous test of SUSY GUT's. 
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10 It should be stressed that this model is quite different from the one discussed in Section 2, since all 
the light Higgs doublets reside solely in the spinor representations 16h + 16h- 
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Figure 3: The proton lifetime (in years) in the dominant channels as a function of £ for 
Strategies I and II. See the caption of Figure ([[]) and Sections fO| and || for the explanation. 



17 



References 

[1] Albright C.H., Barr S.M., Phys. Rev. D62, (2000) 093008 

[2] Aoki S. et al, Phys. Rev. D62, (2000) 014506 

[3] Arafune J., Nihei T., Progr. Teor. Phys. 93 (1995) 665 

[4] Babu K.S., Pati J.C., Wilczek F., Nucl. Phys. B566, (2000) 33 

[5] Buras A.J., Ellis J., Gaillard M.K., Nanopoulos D.V., Nucl. Phys. B135, (1978) 66 

[6] Chadha S., Daniel M., Nucl. Phys. B229, (1983) 105 

[7] Claudson M., Hall L.J., Wise M.B., Nucl. Phys. B195, (1982) 297 

[8] Dermisek R., Man A., Raby S., Phys. Rev. D63, (2001) 035001 

[9] Dimopoulos S., Raby S., Wilczek S., Phys. Lett. B112, (1982) 133 

[10] Dvali G., Pokorski S., Phys. Lett. B379, (1996) 126 

[11] Ellis J., Nanopoulos D.V., Rudaz S., Nucl. Phys. B202, (1982) 43 

[12] Fritzsch H., Minkowski P., Ann. Phys. 93, (1975) 193 

[13] Georgi H., Glashow S.L., Phys. Rev. Lett. 32, (1974) 438 

[14] Goto T., Nihei T., Phys. Rev. D59, (1999) 115009 

[15] Hisano J., Murayama H., Yanagida T., Nucl. Phys. B402, (1993) 46 

[16] Hisano J., |hep-ph/0004266 



[17] Murayama H., Pierce A., |hep-ph/0108104 

[18] Lucas V., Raby S., Phys. Rev. D54, (1996) 2261 

[19] Pati J.C, Salam A., Phys. Rev. D8, (1973) 1240 

[20] Sakai N., Yanagida T., Nucl. Phys. B197, (1982) 533 

[21] Super-Kamiokande collaboration, Phys. Rev. Lett. 83, (1999) 1529 

[22] Weinberg S., Phys. Rev. D26, (1982) 287 



18 



